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Ab initio local-density-functional-theory calculations of 
formation energies, surface stress, and multilayer relaxations 
are reported for the (111), (100), and (110) surfaces of Rh. 
The study is performed using ultrasoft pseudopotentials and 
plane waves in a parallel implementation. 



I. INTRODUCTION 

In this paper we report ab initio calculations of surface 
energies, stress and relaxations of low-index Rh surfaces. 
Our results provide, among others, original information 
on surface stress and surface stress anisotropy, confirming 
the recent findings by Feibelman for Pd and Ft This 
is a positive test of the feasibility and accuracy of such 
calculations within the pseudopotential method. 

A proper understanding of the basic properties of clean 
metal surfaces is an essential step forward in the study 
of adatom diffusion and of catalysis of basic chemical re- 
actions on metal surfaces, both topics of obvious funda- 
mental and technological interest j^. Ab initio calcula- 
tions on transition metal surfaces can provide adsorption 
geometries, dissociation barriers, and diffusion paths for 
use within empirical methods (e.g. kinetic Montecarlo), 
and deeper insight into the basic governing mechanisms, 
and are therefore of the utmost importance. 

With a view at studying diffusion barriers and 
molecular dissociation B on low-index surfaces of tran- 
sition metals, we have undertaken a series of ab initio 
calculations on the surfaces of Rh. Besides its interest as 
a catalyst, good experimental data exist for seldiffusion 
on Rh, in particular for its (100) face Here we present 
surface formation energies, surface stress, and multilayer 
relaxations of the clean Rh low-index surfaces; work in 
progress on surface vacancy formation, homoadsorption 
and self-diffusion on Rh surfaces will be presented else- 
where. Our calculations were performed within local- 
density- functional-theory using ultrasoft pseudopo- 
tentials 1^] to describe ion-electron interaction, a plane- 
waves basis, and iterative diagonalization algorithms in 
an efficient parallel implementation A newly devel- 
oped variant of the Vanderbilt scheme has been used, 
whereby norm conservation is only released for those 
angular-momentum channels which would otherwise gen- 
erate very hard potentials (the c?-channel in the present 
case) 0. 

For bulk Rh, the following computational parameters 
were found to give converged total energy differences 



and structural properties: plane-wave cutoff of 30 Ryd, 
10 special-point mesh for the fee lattice, Fermi-surface 
smearing of 0.05 Ryd, first-order approximation of Meth- 
fessel and Paxton jl^] for the occupation numbers distri- 
bution function (an up-to-date, lucid discussion of this 
treatment of the metallic state in ab initio calculations 
has been given recently by de Gironcoli |ll| ) . For all cal- 
culations, we used downfolded meshes equivalent to the 
one used for the bulk. The resulting bulk parameters 
are qq = 7.215 bohr and B = 3.079 Mbar. The theo- 
retical lattice constant was used in all calculations. The 
surface areas per atom are 26.03 bohr^, 22.54 bohr^ and 
36.81 bohr^ for the (100), (111) and (110) surfaces, re- 
spectively. Our results, given below in eV/atom, should 
be divided by the above values to switch to eV/bohr-^. 



II. RELAXED SURFACE GEOMETRIES 

The topmost layers of most clean transition-metal sur- 
faces are known to relax inward . The intralayer spac- 
ing between the first two top layers, in particular, is ap- 
preciably reduced with respect to its bulk value. Such a 
relaxation is accompanied by smaller relaxations of the 
second and third layers, that may be directed both in- 
ward or outward. 

The magnitude of the relaxations depends markedly on 
surface orientation. It is common knowledge that larger 
inward relaxations occur for rougher surfaces (i.e. with 
lower atomic density). Our results, reported for the three 
low-index clean surfaces of Rh in Table ||, confirm this 
tendency: the (110) surface shows the largest relaxations, 
the closest-packed (111) the smallest. Also, even the sec- 
ond and third layer in (110) relax noticeably, while for 
(100) and (111) only the top-layers show an appreciable 
shift. 

The results in Table | were obtained with 7-layers slabs, 
whereby all layers have been relaxed to their equilibrium 
position. The relaxation pattern is insensitive to a fur- 
ther increase of the number of layer for the (100) and 
(111) surfaces. For the (110), changes are non-zero but 
marginal: using a 9-layers slab, Adi2 is almost unchanged 
(-10.7 %), and Ad23 and Ad34 are only slightly reduced 
(3.5 % and -1.0 % respectively). We emphasize that re- 
sults in Table Q are variations in interlayer spacing, which 
include relaxation contributions from the layers above 
and below the one considered. For the top layer these 
values are equal to the deviation from the ideal position 
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only when the underlying layers are kept fixed. 

In Table |, our results are also compared with those of 
Methfessel et al. who used the FP-LMTO method 
to calculate the top layer relaxations while keeping all 
other layers fixed. While the relaxation trend as function 
of roughness is the same in both calculations, a sizable 
quantitative difference appears for the (110) surface. To 
clarify this, we repeated our calculations relaxing only 
the first surface layer, obtaining a relaxation of -7.5 %. 
The deviation from our full results is thus due to the ne- 
glect of multilayer relaxation in Ref. |jl^. Indeed, this 
relaxation mode is expected to play a role on the (110) 
face, where first-layer and third-layer atoms are nearest- 
neighbors. Our results show that the assumption (sug- 
gested in Ref. for bcc (100) surfaces) that the top 
interlayer spacing remains unchanged irrespective of the 
movements of the layers below, is inappropriate for the 
fee (110) surface. 

We now briefly recall Pettifor's ||l^ model for the trend 
systematization of the basic properties of transition- 
metal surfaces. In transition metals, according to this 
model, the bulk equilibrium lattice constant results from 
the competition of a negative (attractive) pressure due 
to bonding, localized, and directional d states, and of 
the repulsive pressure of the delocalized sp ones (which 
tend to decrease their kinetic energy). At the surface, sp 
electrons can spill out into vacuum, and the top layer 
is further driven inward by the d electron attraction. 
Within this model, one expects the top-layer relaxation 
to follow a roughly parabolic trend as a function of d oc- 
cupation. (This feature was basically confirmed by the 
DFT-LDA calculations of Methfessel et al. Q, where 
further detailed discussions can be found.) At the end of 
a transition series, the s-d balance is reversed, with the 
s electrons providing cohesion and the d states function- 
ing more and more as a closed shell as one approaches 
the noble metals. The fee structure is typical of this 
section of the transition series. Rh is peculiar in that 
it is the first metal in the Ad transition series found to 
be stable in the fee structure, while still exhibiting typ- 
ical transition metal-like features such as sizable inward 
relaxations, high Fermi level density of states (mostly 
d in nature), and the ensuing efficient and short-ranged 
screening. 

III. SURFACE FORMATION ENERGIES AND 
WORK FUNCTIONS 

Surface roughness plays a key role in determining the 
general trends of formation energy, work function, and 
relaxations. One could pick any of these quantities as a 
quantitative measure of the roughness, since each of them 
exhibits well defined trends as a function of roughness. 
Within the Smoluchowsky model ||l^, a rough surface 
presents a smaller inward-oriented electric dipole moment 
compared to a smoother surface of the same material. 



One thus expects the (110) surface to have the largest 
formation energy (having the highest number of removed 
nearest neighbors) and the smallest work function; on the 
other hand, the close-packed (111) surface should have 
the smallest surface energy and the largest work function. 
Our results for the surface energies a and work functions 
W of the three surfaces, listed in Table ||, are in line with 
these expectations. 

Results are given for both relaxed and unrelaxed sur- 
faces. Surface energies and work functions are only 
marginally affected by relaxations, even for the rough 
(110) surface. The discussion of this feature given in a 
previous work on W applies to Rh as well. Our re- 
sults are seen to compare well with those of Ref. [ p^ . 
Numerical deviations are quite minor, considering the 
difference between the two methods employed. 

A technical point concerning the calculation of sur- 
face energies is in order here. It is generally assumed 
that surface energies can be extracted straightforwardly 
as differences of the energy of a slab mimicking the sur- 
face, and of an appropriate bulk energy. For the present 
case 

a= hm i(£;i^^^-iV£;buik), (1) 

with N the number of layers and i?buik the bulk energy 
per atom. It has been recently proven by Boettger |1^] 
that the surface energy diverges as a function of slab 
thickness if the incremental energy of the slab upon ad- 
dition of an atom differs (however little) from the bulk 
energy i?buik- While this results has not yet been widely 
appreciated, in a recent work it has been shown p8[ | that 
the problem is indeed relevant in practice even for tech- 
nically state-of-the-art calculations. A simple solution to 
the problem was suggested : the bulk energy to be 
used in Eq. [I] is the linear part of the slab total energy 
as a function of N. This choice of £^buik, besides being 
the natural one, leads to fast convergence of a (see Ref. 
jlSf for further details). In our calculations we indeed 
encountered essentially the same situation mentioned in 
Ref. for Ft (100), and computed a by the method 
suggested there. 

IV. SURFACE STRESS 

Surface stress is an important quantity providing in- 
sight into surface structure, and useful trend information 
on a number of processes, such as surface recontructions. 
For example, relief of tensile in-plane stress is believed to 
cause the quasi- hex reconstruction of the (100) surface of 
5d fee transition metals [l9| ] ; a contractive reconstruction 
observed recently on Cu (100) has been attributed to 
the same cause The surface stress is defined as the 
strain derivative of surface energy per surface cell |]l|,0 . 
If all bulk contributions to the stress are zero, the surface 
stress is thus 
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where Tq,^ is the volume- averaged stress tensor [ p2| 



of 

the supercell. As usual the factor of 2 accounts for the 
two surfaces of the slab. A positive t^^^ is a tensile 
stress favoring in-plane contraction of the surface, while 
a negative (compressive) stress favors surface expansion. 

For the stress calculations, we used an energy cutoff 
of 45 Ryd, at which the bulk stress is essentially zero 
^^buik ^ 5 X 10^*^ eV/bohr^), and dispensed therewith 
with bulk corrections to the surface stress If the 

bulk stress were non zero, it should be subtracted out 
much in the same way as in the calculation of the surface 
energy. In passing we verified that the surface energies 
calculated at 45 Ry deviate from those calculated at 30 
Ry by less than 0.01 eV/atom. 

In Table [II we report the surface stress for the (111), 
(100) and (110) surfaces of Rh. For the unrelaxed sur- 
face the planar components of the stress are sizably larger 
than the vertical component. The latter is however non- 
zero, and in fact quite large for the very open (110) sur- 
face. Relaxation reduces appreciably the in-plane compo- 
nents (up to about 50 %) and renders the z component of 
the supercell stress negligible in comparison to the planar 
components, as it should be. 

As expected, the stress is anisotropic for the (110) 
surface. The ratio of the stress components along the 
[iTO] and [001] directions is 1.59 (1.09) for the relaxed 
(unrelaxed) surface. These appear consistent with the 
anisotropy values of 1.47 (1.05) for relaxed (unrelaxed) 
Pd (110), and 2.08 (1.47) for relaxed (unrelaxed) Ft (110) 
recently reported by Feibelman 0] . Both anisotropy and 
absolute values of the (110) stress components are compa- 
rable to those of Pd, but much smaller than for Pt. This 
seems consistent with the absence of missing-row recon- 
structions of clean Rh (110) and Pd (110), and its pres- 
ence for Pt (110). Note also that anisotropy is strongly 
enhanced by relaxation due to a larger decrease of the 
stress along [001] (i.e. transversally to the surface chan- 
nels), again in agreement with Feibelman's results for 
Pt and Pd. Also, the stress values are quite insensitive 
to the number of layers used in its calculation, which is 
consistent with Feibelman's finding that close to 99 % of 
the stress on a Pt (111) surface originates from the top 
surface layer. 

To our knowledge these are the first calculations of 
stress for Rh surfaces. The only other stress calculation 
for Rh we are aware of is that of Ref. ||l^, in which a 
stress of 1.94 eV/atom was obtained for unrelaxed Rh 
(100). This compares well with the value of 2.03 eV/atom 
obtained here, expecially in view of the very different 
computational methods used in the two cases. 



V. SUMMARY 

The present ab initio calculations of surface formation 
energies, work functions, relaxations, and stress for the 
clean low-index surfaces of Rh provide an encouraging 
test of the accuracy of the (ultrasoft)-pseudopotential 
method for transition metal surface studies. It is antici- 
pated that it will soon be possible to study e.g. the sur- 
face stress of reconstructed surfaces, which should clarify 
the role of stress and stress anisotropy in surface recon- 
structions. The study of adsorption and of simple surface 
defects (vacancies, adatom- vacancy pairs), and diffusion 
processes, also seems within reach. 
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TABLE III. Surface stress (eV/atom) for unrclaxod ('u') and 
relaxed ('r') Rh surfaces. For the (110) surface, the xx component 
is the [001] one, and the yy is the [110]. See the text for the effects 
of relaxation. 
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TABLE L Intra-layer relaxations (percentage variation with re- 
spect to ideal layer spacing) for the three low-index surfaces of Rh. 
Results from Ref. [h3] are included for comparison. 
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TABLE H. Surface energy cr, relaxation energy Act (the surface 
energy change upon relaxation) and work function for the three 
clean surfaces, 'u' and 'r' indicate unrelaxed and relaxed surfaces 
respectively. All results are in eV/atom. Results of Ref. [Q are 
included for comparison. 
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